Entropic Splitter for Particle Separation 
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We present a particle separation mechanism which induces motion of particles of different sizes 
in opposite directions. The mechanism is based on the combined action of a driving force and an 
entropic rectification of the Brownian fluctuations caused by the asymmetric form of the channel 
along which particles proceed. The entropic splitting effect shown could be controlled upon variation 
of the geometrical parameters of the channel and could be implemented in narrow channels and 
microfluidic devices. 
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Matter very often manifests as divided into very fine 
parts whose nature dictate the overall properties of the 
system. Grains, conglomerates or meso-structures dis- 
persed in a liquid or liquid-like phase are basic structures 
frequently found in many physico-chemical and biological 
problems. The dissimilarity of these units, resulting from 
random fractionation or self-assembling processes giving 
rise to a very disparate distribution of sizes, originates 
a heterogeneous response of the system which makes its 
typification difficult. To obtain pure substances by sep- 
arating wanted from unwanted elements presents a main 
challenge in basic research and industrial processing, and 
in nanotechnology as well. 

Particle separation techniques use the fact that the re- 
sponse of the particles to external stimulus, such as gra- 
dients or fields, depends on their size. Filtering particles 
of different size is traditionally performed by means of 
centrifugal fractionation [lj, phoretic forces [2|t4| or ex- 
ternal fields [5j- By means of these methods, the sorting 
of particles proceeds either by size exclusion, as happens 
in a sieve, or by migration through the host medium, a gel 
or porous media. In these cases all particles move in the 
same direction but at different speeds. Novel separation 
techniques based on flashin g B -|a, drift "determinis- 
tic" jTo| and geometric [T3jl2| ratchets have also been 
proposed for sorting [l3| . 

Here we present a novel splitting mechanism that sep- 
arates particles in different directions by purely entropic 
means. The working principle relies on the combined ac- 
tion of a static force and an entropic rectification fl4- 17|. 
Small particles follow the force whereas the motion of big 
particles is rectified to proceed in the opposite direction 
resulting in a faster splitting. This mechanism could be 
readily implemented in microchannels or microfluidic sys- 
tems. The geometry of the channel can be tuned to be 
very selective leading to an efficient separation of parti- 
cles of very similar radius. 

The Brownian motion of particles in confined geome- 



tries exhibits a very rich and striking phenomenology 



23j . The effects of confinement can be described by 
means of an effective entropic potential resulting from 
the variation of the space accessible to particles along 
the transport direction 18, 24, 25)]. The height of the en- 
tropic barrier associated to the bottlenecks depends on 
particle radius. In an asymmetric channel, the resulting 
entropic potential becomes asymmetric, and rectification 
of a zero-mean oscillating force can occur. The strength 
of this rectification depends on the particle radius and is 
stronger for large particles (for which the entropic barrier 
is larger, see Fig. [T|). Thus, in the presence of a small 
static force directed in the opposite direction of recti- 
fication, it is possible to separate particles of different 
sizes. The trajectories of small particles mainly follow 
the force, whereas large particles move in the opposite 
direction. To illustrate this effect we have chosen the 
geometry depicted in Fig. [T] 

The dynamics of a Brownian particle in a 2D channel 
as the one depicted in Fig. [TJ driven by a static force / 
and an oscillating (square wave) force F(t), both applied 
along the principal axis of the channel, can be described 
by means of the Langevin equation which, in the over- 
damped limit, reads 



dr 



7 ^ = - (/ + F(t)) ex + Vlk B TZ(t) , 



(1) 



where r denotes the position of the particle, 7 is the 
friction coefficient, e x is the unit vector along x-direction, 
and is a Gaussian white noise with zero mean which 
obeys the fluctuation-dissipation relation £j(t')) = 
2 5ij5(t — t') for i,j — x,y. The explicit form of the 
oscillating driving force is F(t) = A sgn [sin(fii)] where A 
is the amplitude, sgn[£] represents the sign function, and 
O is the driving frequency. This choice does not represent 
a restriction on the particle splitting effect which could 
also be observed for a sinusoidal driving force. 

The Langevin equation (JXJ) must be solved by imposing 
vanishing outflow at the walls of the structure. For the 
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FIG. 1: (top) Schematic illustration of the two-dimensional 
channel confining the motion of the Brownian particles. The 
structure is defined by Eq. ([2]). Brownian particles are driven 
by a constant force / and a square wave force F(t) along the 
longitudinal direction. The values of the chosen parameters 
defined in Eq. © are: m 1 = 2, m 2 = 0.4, b/L = 0.1. The 
dashed lines represent the limit for the positions of the cen- 
ter of the particles within the channel, (bottom) Effective 
entropic potential for the two particle sizes depicted above. 



2D structure depicted in Fig. [TJ the walls are denned by 



Vu(x) 



m 2 (L 



if x < c 
otherwise 



(2) 



where yu(%) and yi(x) = —y u (x) correspond to the 
upper and lower boundary functions, respectively, b is 
the half-width of the bottleneck, mi and mi are the 
slopes of the walls, L is the periodicity of the channel, 
c = Lraij(m\ + ra-i) indicates the location of the point 
of maximum width, and x = x mod L is the modulo 
function (to create a periodic structure) cf. Fig. [TJ 

For a hard particle of radius r inside the channel, the 
space available for its center is restricted by a distance r 
from the walls, reading 



w u (x) = 



~\Jr 2 — x 2 + b, 
b + m\x — ry (1 + mf), 
b + m 2 {L-x) -r yfiT 
-y/r 2 -(x-L) 2 + b, 



< x < o p 

Of) X Ct) 



>>2J> 



(3) 



< x < L 



where o p = rmi/y (1 + m 2 ), L p = L — rm^/ y(T+ mf), 
and c p = c + r (^J (1 + m\) - /( TO i + m ^)- 

The lower parallel curve is just w\(x) = —w n (x). Conse- 
quently, 2w(x) — w u (x) — wi(x) gives the local width of 
the structure accessible for the center of a hard particle 
of radius r. Assuming throughout a dilute particle den- 
sity and a strong viscous, low Reynold number dynamics, 
all relevant hydrodynamic wall-particle interactions are 
small and of repulsive character [26j. Effectively this re- 
sults in a slightly increased effective particle radius; i.e., 



r e ft > r. Such hydrodynamic interactions thus result in 
somewhat larger average transport velocities, note Figs. 
2a, b below. 

Also the friction depends on size, and using Stokes' 
law as an approximation one obtains 7 = "for/b and 
D = D^b/r, where 70 and Dq = k^T/^o are the friction 
and diffusion coefficients of a particle of radius equal to 
the bottleneck half-width b. To mimic the particular case 
of DNA electrophoreses, we will consider that the forces 
depend linearly on the radius of the particles. Specifi- 
cally, we set 



/ = f r/b and A = F r/b , 



(4) 



where /o and Fq are the strengths of the static and peri- 
odic forces for a particle of radius r. 

For the sake of a dimensionless description, we scale 
all variables using three characteristic parameters: the 
characteristic length L, energy k^T, and diffusion coef- 
ficient Dq. Particularly, x = x/ L, y = y/L, b = b/L, 
w\ = Wi/L = —w u ; t = tjr and 51 = fir, where 
r = L 2 / Dq is the characteristic diffusion time. The 
scaled forces are: / = fL/k B T and F(t) = F(t)L/k B T. 
In the following we shall omit the tilde symbols. In di- 
mensionless form the Langevin equation ([l} reads: 

dr 
dt 



= - (/o + F sgn[sin(Qt)]) e x + yjb/r&t) . (5) 



The description of this system alternatively can be 
given by the concept of an entropic potential 11 81 1 and the 
corresponding Fick- Jacobs (FJ) equation [IF 




et dx\ D ^{^ +v ^ p 



(G) 



where 



V(x) = U- TS(x) = [F(t) +f]x- In [2 w{x)} , (7) 

is the free energy including an entropic contribution 
TS(x) = In[2to(as)], D(x) = b/{r[l + w'ix) 2 } 1 / 3 } is the 
position-dependent diffusion coefficient, and the prime 
refers to the derivative with respect to x. This approxi- 
mation is expected to be very accurate for bias strengths 
|/o I < 1) amplitudes F < 1 and for w'(x) 2 < 1 |21|. 

In the adiabatic limit, the average velocity can be cal- 
culated as 



J(F ) + J(-Fq) 



(8) 



where the current J(Fq) is given by Stratonovich's for- 
mula (UEl. 



J(F ) = 



1 



-(F0+/0W& 



(9) 



Note that for very large driving strengths Fq 3> 1, while 
keeping a finite bias strength |/o|, i.e. far beyond the 
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FIG. 2: Average current vs. the amplitude of the periodic 
forcing Fo in the adiabatic limit ( Q = flo = 7r/10) for parti- 
cles of different radius, for a channel with mi = 2, m2 = 0.4. 
The solid lines indicate the results of the simulations, whereas 
the dashed line plots the prediction according to the FJ equa- 
tion for r = 0.9b. For /o < small particles go to the left, 
whereas big particles move to the right. The critical size that 
divides these two behaviors can be tuned by the value of Fq. 



regime of validity of the Fick-Jacobs approximation, the 
average velocity approaches the deterministic behavior 
lBL 19| . In this limit we then find that J(F ) ~ — J(— Fq) 
so that an asymptotic vanishing net velocity < v > 
emerges. Note that this asymptotic regime of very large 
amplitude strength Fo is not yet reached in Fig. 

Fig. [2] plots the average current vs. the amplitude of 
the periodic forcing Fq, in the adiabatic limit (Q = 7r/10), 
for particles of different radius. In absence of static forc- 
ing, this asymmetric channel rectifies the oscillating force 
giving rise to a net positive velocity, i.e. all particles 
move towards the right (see Fig. [5Ji). Its magnitude de- 
pends on the strength of the rectification which is more 
intense the larger is the particle radius. If we apply a 
small static force in the negative direction (see Fig. 
for intermediate periodic forcing strengths Fo, particles 
larger than a given threshold radius move to the right, 
whereas particles smaller than that move to the left. In 
this way, one can separate particles of different radii and 
make them move in opposite directions. The splitting 
effect is illustrated schematically with Fig. [3] 

Fig. [5] also represents the results for the average veloc- 
ity obtained from the FJ equation for the case r = 0.96. 
The agreement with the simulation results is very good at 
small values of the force (for Fq < 10). At larger values of 




FIG. 3: Schematic illustration of the functioning of the en- 
tropic splitter. A mixture of particles of two different radii 
is placed initially in the center. Under the combined pres- 
ence of the static and periodic forcing, large particles move 
towards the right, whereas small particles follow the bias, i.e. 
are moving towards the left. 








FIG. 4: Average current vs. the value of the static force fo 
in the adiabatic limit (Q — n/10) for particles of different 
radius, for the channel plotted in Fig. [ljwith Fo = 20. In all 
cases the velocity increases linearly with fo. As fo gets more 
negative, the particles invert their velocity and start to move 
to the left. The critical force that determines this velocity 
inversion depends on the radius r and gets progressively larger 
(in absolute value) for larger radii. 



Fo, the agreement is better when the radius increases or 
equivalently when entropic effects are more pronounced. 
In the domain of validity of the FJ description the cor- 
responding applied forces already yield a significant en- 
tropic splitting. 

The splitting effect can be tuned by the value of ei- 
ther the amplitude of the periodic forcing or the static 
force. Fig. [4] represents the average value of the veloc- 
ity vs. fa for the same geometry and Fo = 20. One 
can see that by tuning fo, one can control the separation 
of particles of different sizes. For instance, by selecting 
fo = —1.5, small particles of radius r = 0.16 will move to 
the left with a velocity —0.7, whereas large particles of 
radius r — 0.96 will drift to the right with velocity 0.5. 
More importantly, the velocity depends almost linearly 
on fo, thus facilitating an efficient control of the sepa- 
ration effect. In addition, by progressively changing fo, 
one obtains a device that, with a fixed geometry, can be 
used to continuously separate particles of any size. 

We have also analyzed the effect of the frequency of 
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the periodic forcing on the particles current. For small 
frequencies the results agree with those obtained in the 
adiabatic limit. As the frequency increases the velocity 
becomes progressively smaller. Eventually, at very high 
frequencies, the change in the oscillating force is so fast 
that the particle cannot follow it and as a result there is 
a vanishing effect of the oscillating force. 

The fact that particles of different sizes travel in op- 
posite directions leads to an efficient separation that can 
be improved further by increasing the number of peri- 
odic cavities. A way to quantify this is to calculate the 
probability that a large particle with an average positive 
velocity v > reaches the "wrong" collector, i.e. the 
left boundary of the device placed at x — nL, where n 
is a negative integer (see Fig. [3]). This probability is 
P(x < nL,t) = 0.5 + 0.5 erf ((nL - vt)/y/4D eS t) where 
we have assumed that the process can be described as a 
driven-diffusion process with an average velocity v and an 
effective diffusion coefficient D c g. This probability peaks 
in time at t max = —n/v, and its maximum value is at 

P(x < nL,t max ) — 0.5 + 0.5 erf f— \/v\n\/D e &\. We use 
this maximum probability as a measure of the purity of 
the sample collected at the end of the device. For typical 
values of the velocity v ~ 0.2 (cf. Fig. 2]) and D e ff < 1, 
one can achieve 99.9968% of purity after only 40 periods. 

Note, that the entropic splitting effect can be used for 
separation of DNA fragments of different size using typ- 
ical values as for DNA electrophoresis. The DNA frag- 
ments are considered as randomly coiled polymers with 
radius R = d(M /Mq) 1 ' 2 where M is the number of base- 
pairs, d = lOOnm is the Kuhn length, Mg = 300 is the 



number of base pairs per Kuhn length [llj . Assuming 
Stokes' law and the Sutherland-Einstein relation, (i) the 
friction coefficient for the dynamics of such a DNA frag- 
ment in water is rj « 2 ■ 10~ 9 kg/s \J M/Mq and (ii) the 
diffusion coefficient is D ■ 
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2 / Sx /Mo7m. 



For a device with L = 3/j,m and the same geometry 
and parameters as in Fig[TJ typical values of the velocities 
are L/t = Do/L ~ 0.9/im/s, achievable by electrophore- 
sis with an electric field of 7 V/cm. For DNA chains of 
R ~ 250nm differing in radius by a ~ 25%, a purity of 
separation of 99.997% results after 40 periods. These val- 
ues compare favorably with commonly used techniques. 
Much larger velocities and efficiencies are possible with 
smaller DNA chains, or nanosized particles. 

The effectiveness of the entropic splitting can be tuned 
by choosing the geometry of the channel. In particu- 
lar, the entropic splitting effect becomes more important 
upon increasing the asymmetry of the walls, the slopes 
of the channel or by decreasing the bottleneck width b. 
Altering the design of the channel geometry it is feasible 
to separate particles of very similar radii, e.g. by choos- 
ing b close to r. Moreover, yet another advantage of this 
set up is that it works in a time-continuous mode. 

In summary, we have presented a novel, purely entropic 



particle splitting mechanism which is able to separate 
particles of different sizes. The mechanism is based on 
the presence of an entropic rectification of fluctuations 
caused by the asymmetric form of the channel. This rec- 
tification may overcome the effect of an applied force by 
reversing the motion of the particles. The mechanism is 
very efficient and can be controlled by tuning the geomet- 
ric parameters of the channel leading to both, different 
velocities and directions. This idea could be implemented 
in constrained structures with narrow channels and pores 
where entropic effects are important. 
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